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ABSTRACT 

To understand the origin of Solar Energetic Particles (SEPs), we must study their injection time 
relative to other solar eruption manifestations. Traditionally the injection time is determined using 
the Velocity Dispersion Analysis (VDA) where a linear fit of the observed event onset times at 1 AU to 
the inverse velocities of SEPs is used to derive the injection time and path length of the first-arriving 
particles. VDA does not, however, take into account that the particles that produce a statistically 
observable onset at I AU have scattered in the interplanetary space. We use Monte Carlo test particle 
simulations of energetic protons to study the effect of particle scattering on the observable SEP event 
onset above pre-event background, and consequently on VDA results. We find that the VDA results 
are sensitive to the properties of the pre-event and event particle spectra as well as SEP injection 
and scattering parameters. In particular, a VDA-obtained path length that is close to the nominal 
Parker spiral length does not imply that the VDA injection time is correct. We study the delay to 
the observed onset caused by scattering of the particles and derive a simple estimate for the delay 
time by using the rate of intensity increase at the SEP onset as a parameter. We apply the correction 
to a magnetically well-connected SEP event of June 10 2000, and show it to improve both the path 
length and injection time estimates, while also increasing the error limits to better reflect the inherent 
uncertainties of VDA. 

Subject headings: Journal approved keywords 


1. INTRODUCTION 

During solar eruptions, charged particles are acceler¬ 
ated up to relativistic energies, to form the solar ener¬ 
getic particle (SEP) population of the cosmic rays ob¬ 
served by in situ instruments at different locations in the 
heliosphere. The particles are believed t o be accelerate d 
in flares and CME-driven shock waves (Reames 1999|). 
However, the relative importance of the hare and CIVIE 
processes on the origin of the observed SEP populations 
is still under scientific discussion, and opinions differ on 
how the eruption phenomena and the SEP production 
are connected (e.g.,|Can e et al.||2010[ |Gopalswamy et al. 
2012 Aschwandenf 2012 1) . 

The difficulty in deducing the particle acceleration sce¬ 
narios during solar eruptions stems from the nature of 
SEP observations. The propagation of the charged SEPs 
is affected by the interplanetary magnetic field. The 
SEPs are guided by the large-scale Archimedean spiral 
structure of the interplanetary magnetic field, the Parker 
Spiral. Solar wind is turbulent, and the particles scatter 
off th e inhomogeneities of the magnetic field (e.g., Parker 
19651. Thus, the particle propagation is diffusive rather 
than direct propagation from the acceleration site to the 
in-situ particle detectors, typically at 1 AU from the Sun. 
In order to understand the connections between the com- 

timo.lm.laitinen@gmail.com 


ponents of a solar eruption and the observed SEP intensi¬ 
ties, we must understand the propagation of SEPs in the 
interplanetary space, and deconvolve it from the obser- 
vati ons. This has been done for several SEP events (see, 
e.g., Kallenrode||I993[ ITorsti et al.||1996| Laitinen et al. 


2000 


2003 Agueda et al.|2069p , with most recent 


works introducing c ross-held diffusion into the modelling 
of SEP events (e.g.,|Zhang et al.|2009[[Droge et al.|2010 



Deconvolving of the interplanetary transport from the 
in-situ SEP observations is, however, not simple, and is 
usually performed only in case studies. For larger sta¬ 
tistical studies, simpler methods to obtain the injection 
time of the SEPs are commonly used. The popular choice 
is to use the Velocity Dispersion Analysis (VDA), where 
the first-observed particles are assumed to have propa¬ 
gated without scattering (e. g. Lin et al.|[l981[ |Reames 


et al. 19851 Torsti et al.| 19981 iKrucker & Lin|2ubu| I'l'ylka 

et al. 2003 Reames||2009 Vainio et al. 2013p. The pos¬ 

sible uncertainties on the arrival times ot hrst particles 
caused by interplanetary transport effects are typically 
not evaluated in these studies. 

The validity of the VDA method has been studied by 
using numerical SEP simulations that solve the focused 
transport equation describing particle propagation in in¬ 
terplanetary space. By using onset times at 1 AU ob¬ 
tained from simulated SEP time-intensity profiles at dif- 
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ferent energies, it has been shown that for strong scat¬ 
tering conditions the VDA can result in large errors 


with the particle distribution function, n, /r, 

z and t the speed, pitch angle cosine, position along 


cles (iKallenrode & Wibberenz 

1990 

Lintunen & Vainio 

V the scattering frequency, and Q describing the parti- 

2UU4 

yaiz et al. 

2DDS 

). i'hese studies, however, dehned 

cle source. The equation is solved using Monte Garlo 


the time when the intensity reaches 1% of the maxi¬ 
mum intensity), which is not the common practice when 
analysing real SEP events. Unlike in simulated SEP 
events, real events may have a complex structure, due to 
local interplanetary magnetic field structures and multi¬ 


ple, energy-dep endent injection components (e.g. Laiti¬ 


nen et al. 


2000). This may affect the time profile of the 


event before the maximum intensity is reached. Thus, 
to get the best estimate of the beginning of the injec¬ 
tion, the onset of an SEP event is typically determined 
by using the moment when the SEP intensity exceeds 
the pre-even t background by a statistically si gnificant 


et al. 


amount (e.g ., Huttunen-Heikinmaa et al. 2005). 


Laitinen 


(2010) analysed simulated SEP events with V UA 


by adding a pre-event background to the time-intensity 
evolution of the SEP event, and using the background- 
exceeding time as the event onset time. They found that 
the errors in VDA have strong dependence on the pre¬ 
event and event maximum spectra. 


In this work, we build on the results of Laitinen et al. 


(2010) and discuss the complexity of the effect of the 


pre-event background and the role of scattering in differ¬ 
ent types of SEP events, and show how different event 
types can result in very different VDA results. We study 
the SEP time-intensity profile of energetic protons at the 
time of the observed onset, and derive a simple estima¬ 
tion for the delay the particles experience due to the 
scattering in interplanetary space. We show that this 
delay estimate can successfully be used as a correction 
to the observed onset times to improve the accuracy of 
the VDA. We apply the correction and its error limits 
to simulated SEP events and the SEP event of June 10 
2000, and show that the method improves the deduced 
injection time in both cases. 

To clarify the used terminology, we use the term ‘‘‘'in¬ 
jection time" to describe the time of release of SEPs at 
or near the Sun, “onset time" to describe the time when 
the SEP intensities are observed to rise at 1 AU (see dis¬ 
cussion below), and “launch time" as the time when a 
CME is estimated to lift off. Furthermore, the symbol 
ti refers to the injection time obtained by VDA, and to, 
the onset time obtained from the simulated events. The 
time required for a scatter-free particle with velocity v to 
propagate distance s, i.e., s/v, is referred to as “scatter- 
free time". 

2. MODELS 

2.1. Solar particle event modelling 

In this work, we study the effect of pre-event back¬ 
ground on SEP injection times as determined by the 
VDA method, in the presence of interplanetary scatter¬ 
ing and different injection profiles. We solve the focused 
transport equation of energetic protons, 

df d f, 9 ^ 2\ f 


cused along the Parker spiral in a frame co-rotating with 
the Sun, and scattered isotropically in pitch angle, with 
ly = v/(2X ), where A is the parallel scattering mea n free 
path (e.g., Torsti et al.|1996 Kocharov et al.|19^ ). Adi¬ 


abatic deceleration is taken into account by scattering 
the particles in the co-rotating solar wind frame. Propa¬ 
gation of paricles across the Parker spiral direction is not 
considered. 

We use the initial energy range from 1 to 120 MeV with 
a power law spectrum oc E~^. This power law is used 
for all runs, to reduce the parameter space. The particles 
are injected at 0.01 AU, at time t = 0, with a reflecting 
boundary at the Sun, and followed for 48 hours. The 
Parker spiral is parametrised by a constant solar wind 
velocity of 400 km/s and solar rotation period of 25.35 
days, resulting in Parker spiral length of 1.17 AU. We use 
constant radial mean free path, with values = 0.3 AU 
at IGV rigidity to represent moderate scattering condi¬ 
tions, and 1 AU at 1 GV rigidity to represent weak scat¬ 


tering conditions (e.g. Palmer[|1982 ). The mean free path 
is taken to depend on the particle rigidity, as con¬ 

sistent witliQuasilmear theory for Kolmog orov slab spec - 
trum ( Jokipii|[i966 ) and observations (e.g. Droge 2000). 


The Monte Garlo simulations are used to obtain a re¬ 
sponse, as observed at 1 AU, for an impulsive injection 
of protons, at 20 logarithmically spaced energy channels 
between 1 and 100 MeV. The upper limit of the channels 
is chosen to be 100 MeV rather than 120 MeV, because 
particles experience adiabatic deceleration. In order to 
mimic more realistic particle release scenarios, we con¬ 
volve the impulsive responses with two injection profiles 
at the Sun: a fast injection model, with 30 minutes of 
linear increase of injection strength, followed by 570 min¬ 
utes of linear decay; and a slow injection model, with 570 
minutes of linear increase, followed with 30 minutes of 
linear decay. These profiles were selected because CME- 
related SEP acceleration has been estimated to be most 


efficient when the CME is at 5-15 R© (e.g. Kahler 


1994). 


A CME with velocity 2000 km/s reaches 5 Kq m halt an 
hour, corresponding to our fast injection model. The 
slow injection model accounts for slower CMEs, higher 
maximum injection heights, poor connection between the 
spacecraft and acceleration region and possible cross-field 
transport of SEPs, that may cause a gradual injection of 
particles. 

We define our study parameters using SEP intensities 
that can be observed directly at 1 AU. We scale the con¬ 
volved time-intensity profiles of simulated SEPs at 1 AU 
so that the spectrum calculated from the maxima of the 
scaled intensities at each channel forms a power law. 




OeR 


E 


( 2 ) 


where Rev is the SEP event maximum intensity at en¬ 
ergy Eq = 88 MeV, and aev the power law index of the 
spectrum. The resulting event time-intensity profiles are 
overlaid on a constant pre-event background, which also 
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Fig. 1.— An example of the event maximum spectrum (red dot¬ 
ted curve),the pre-event background spectrum (green dash-dotted 
curve), with the total intensity spectrum at event maximum given 
by solid black curve. The dashed vertical black line depicts the 
energy Ecut-, where the event maximum intensity is an order of 
magnitude times the background intensity. The event maximum 
spectrum power law index is 3, and the background power law 
index 1. 

follows a power law energy spectrum, 

/ E\ ““‘s 

hg = lobg 

where lo^g is the background intensity at Eq and abg the 
background power law index. An example of the spectra 
is shown in Fig. 

2.2. Velocity dispersion analysis 

The Velocity Dispersion Analysis is based on the as¬ 
sumptions that the first particles observed at a given 
distance from the Sun have been released simultaneously, 
propagate the same path length, and experience no scat¬ 
tering or energy changes. Under these conditions, the 
arrival time, to, of the particles to the observer at dis¬ 
tance s along the magnetic field line is given by 

toiv) = U + -, (4) 

V 

where ti is the particles’ injection time at the Sun, s 
the travelled distance and v the particle velocity. Thus, 
knowing the observed onset times at 1 AU, and the ve¬ 
locities of the particles, a simple linear fitting of this 
data according to the Eq. 0 gives the particles’ injec¬ 
tion time at the Sun, and the path length traveled by 
the particles. In the simulated events, the injection of 
the particles starts at t = 0, thus a successful VDA fit 
would give ti = 0 and s = 1.17 AU. 

While all of the VDA assumptions can be questioned, 
in this study we concentrate only on the effect of the scat¬ 
tering on the derived injection time and path length. In 
particular, we will study the common practice of defin¬ 
ing SEP onset as the moment the intensity is discernible 
from the pre-event background, and how this practice 
affects the VDA results. 

The SEP event onset time, observed by particle detec¬ 
tors, is often difficult to determine due to low counting 
rate of particles of the ambient energetic particle popu¬ 
lation before and at the very beginning of the event. To 



Fig. 2.— An example of a proton event with energies 1—100 MeV 
overlaid on a pre-event background. The lowest-energy intensities 
are on the top of the figure, the highest on the bottom. The onset 
times, as defined in the text, are marked by squares. The red 
curves and the large red squares represent the channels where the 
event maximum intensity is over 10 times the bacground, the green 
curves and the small green squares represent the smaller-intensity 
energies. The thick brown curve represents the scatter-free arrival 
time of particles. 



Fig. 3.— An example of the linear fit of the observed onset times 
and the inverse velocities. The squares depict the observed onset 
times, with the color scheme the same as in Fig.[^ The thin black 
curve depicts the resulting linear fit to the energy channels with 
event maximum intensity over 10 times the pre-event background. 
The obtained parameters of this fit are given on the top-left corner 
of the figure. The dashed blue line depicts the velocity dispersion 
for particles that are observed after time t = 1.17 AU/u from their 
injection. 


determine the onset time, typically a threshold of one 
or several standard deviations above the pre-event back¬ 
ground is used to d efine the onset in a statistica.Uy sig nif- 
icant way (see also Huttunen-Heikinmaa et al.||2005 for 
an alternative method), 'i'his results in a delay tor the 
observed onset. In order to mimic this effect, we define 
the onset to be observed when the intensity at 1 AU rises 
10% above the pre-event background. 


3. RESULTS 
3.1. Parametric study 

We show an example of an analysed event in Eigs. 
and In Fig. we show the time evolution of a sim- 
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ulated event, with maximum intensity spectrum cx E~^, 
taking place when the pre-event background spectrum is 
proportional to with the event intensity an order 

of magnitude below the background intensity at energy 
88 MeV, and A,. = 0.3 AU, and fast injection. The onset 
times of the intensity increases are shown by symbols, 
with the brown thick curve crossing the time-intensity 
curves on the left showing the theoretical time of arrival 
of particles in the case of particles following the assumed 
1.17 AU Parker Spiral field without scattering. As can 
be seen, the intensity increase begins significantly later 
than the scatter-free time at all energy channels. 

The effect of the delay on the VDA can be seen in 
Fig. where we show the onset time plotted against 
l/P ^he red and green symbols). As can be seen, the 
onset times are clearly delayed from the velocity disper¬ 
sion pattern for particles that are observed immediately 
after scatter-free propagation to 1.17 AU (dashed blue 
line). In addition, not all observed onsets follow the lin¬ 
ear velocity dispersion pattern. In our example event, 
at large energies the arrival time dependence on 1/j3 is 
inverted, with onsets at higher energeies observed later 
than at lower energies, due to the pre-event background. 
Such energy channels clearly must be excluded from the 
velocity dispersion fit. In order to achieve this, we have 
excluded energy channels where Iev{E)/I^{E) < 10 (the 
green curves and small symbols in Figs. [^andl^. 

However, even after removing the nonTinearly behav¬ 
ing energy channels from the velocity dispersion fit, the 
pre-event background and scattering still have a substan¬ 
tial effect on the determination of the SEP injection time. 
The observed onset time at different energy channels is 
delayed by 20-40 minutes from the scatter-free arrival 
time. The VDA-fitted injection time for this event is 21 
minutes later than the actual release time of particles at 
the Sun. The path length, 1.26 AU, is clearly longer than 
the nominal 1.17 AU Parker spiral length in our model. 
The increased path length is due to the low-energy (high 
l/P) particles having longer delay relative to the scatter- 
free arrival time (blue dashed line in Fig.|^, as compared 
to the higher-energy particles. 

In Fig. we show the VDA fits for A^ = 0.3 AU for 
the fast (left column) and slow (right column) injection 
model, and different types of events. In the top pan¬ 
els, we show an SEP event with the maximum spectrum 
softer than pre-event background spectrum. This type of 
event may be observed when the pre-even t background is 
not a ffected by preceding SEP events (e.g. 


Valtonen et al. 


20011. In such an event, the intensities at higher energies 


are masked more efficiently by the pre-event background 
than at lower energies. As can be seen, the VDA path 
lengths are typically close to or slightly longer than the 
nominal Parker spiral length, with the obtained injection 
time considerably later than the actual solar injection 
time. 

In the event shown in the middle panels, the pre-event 
and event maximum spectral indices are equal. This can 
take place for example when a soft-spectrum SEP event 
follows a harder-spectrum SEP event. The delay with re¬ 
spect to the scatter-free arrival time is signihcantly larger 
at low energies. This results in a significantly longer path 
length, of 1.52 AU in the case of fast injection model 
(left middle panel), and 1.75 AU for the slow injection 
model (right middle panel). The long path length, how¬ 


ever, compensates partly the delay caused by scattering, 
resulting in only 2-minute error for the injection time 
determination with the fast i^ection model. 

The bottom panels of Pig. depict an event where a 
hard SEP event takes place on a soft pre-event back¬ 
ground. This typically takes place when a previous SEP 
event is still decaying at the time of a new SEP event. 
In this case, the lower energies are efficiently masked by 
the pre-event background, and the velocity dispersion fit 
results in a very long path length. As can be seen in 
both the fast and slow injection models (bottom left and 
right panels, respectively), the obtained injection time 
can precede the real injection time in such cases. 

Fig . I^shows that also the injection profile of the parti¬ 
cles has a significant effect on the VDA results, with the 
comparison of fast and slow injection prohles on the left 
and right columns, respectively. As can be seen, the ef¬ 
fect of the pre-event and event maximum spectral shapes 
on the velocity dispersion pattern is similar for both in¬ 
jection models. However, the delay times for individual 
energy channels are longer in the slow injection case, and 
the resulting error in injection time is also larger. 

As shown in Pig.[^ the assumption of scatter-free prop¬ 
agation can result in significant systematic errors for the 
SEP injection time and the traversed path length. This 
error is caused by the delay of particles due to interplan¬ 
etary scattering, which should be taken into account in 
the velocity dispersion analysis as 

g 

toiv)=U-\ - \-tdiv,X,to), (5) 

V 

The estimation of the delay time td is not trivial, as it 
depends on scattering conditions, injection profile and 
the level of the pre-event background relative to the SEP 
event intensity. Also the estimation of the resulting error 
in VDA is not straigthforward: as shown in the middle 
panels of Fig. an error in the path length can compen¬ 
sate the error in the injection time determination. Such 
a complicated relation between the observed velocity dis¬ 
persion and the injection is difficult to analyse, and may 
skew the results of large statistical SEP studies, where 
fitting of the transport of SEPs is not feasible. 


3.2. Analysis of the delay time in VDA 

The rise time of SEP intensities above the pre-event 
background is not the only observable of an SEP event 
onset. As shown in Fig. the intensities rise initially 
roughly exponentially, with different rise rates at differ¬ 
ent energies and intensities. This gradual rise, as opposed 
to immediate rise to the maximum intensity, is a result 
of the particle scattering in interplanetary space. In the 
following, we study the rate of intensity increase at the 
time of SEP event onset. 

Eor a diffusion process with scattering mean free path 
A, the intensity of particles with velocity v at distance r 
and time t is given as 

1 f 3r^ j 

{4Xvt/if^ I 4AutJ 

From this, we can obtain the timescale of the intensity 
increase as 


1 d/diff _ 3 1 
Idiff dt 2 t \ 2 Xvt 


( 7 ) 
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Fig. 4.— Onset times fitted for simulated events with different event maximum and background spectra, and injection profiles. The left 
column presents the results of the fast injection model, and the right the slow injection model, with Ar = 0.3 AU. The symbols represent 
the obtained onset times at different energies, with the large [small] symbols for channels where the intensity at event maximum is more 
[less] than an order of magnitude above pre-event background. The solid curve shows the l//3-fit to the onsets, and the dashed curve the fit 
expected for scatter-free particles. The inset in the top left of the panels shows the event maximum spectrum (dashed curve), background 
spectrum (dot-dashed curve) and the total spectrum (solid curve). 
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Fast onset, \ =0.5 AU 
Slow onset, X'=0.3 AU X 
Fast onset, =1 AU O 

Slow onset, X' = 1 AU 
_ _ _ c _ 
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Fig. 5.— Scatter plot of = ti — s/v against the scaled intensity 
increase scale time, for simulations with different mean free paths 
and inje ction profiles. The solid curve represents the best fit of the 
Eq. ( |10[ | to the simulation data, with the dashed curves depicting 
the error limits of the fit. 

The diffusion approach is valid only for times t ^ A/3z; 
and t ^ r/v. Indeed, the intensity maximum takes place 
at tmax = f/( 2A) r/u, which for A > r/2 is before the 
particles with velocity v can arrive scatter-free to dis¬ 
tance r, clearly an unphysical result. However, as there 
are no analytic descriptions for the initial phase for dif¬ 
fusively spreading particles, we will use diffusion as a 
starting point of our study. 

We are only interested in the initial increase of the 
SEP event, and use only times when the intensity is 
less than /e„(i?)/10. Thus we only use the first term 
in Eq. Q. Rewriting r/ using the diffusion timescale, 
td = 3r^(2Ar;), we arrive with a simple scaling between 
the time from SEP injection at the sun, t, and the inten¬ 
sity increase timescale r/ as 


1 


-titd = t . 


( 8 ) 


To study this scaling, we simulated a large number of 
SEP events with the ratio loev/Iobg ranging from 10“^ to 
10^, and the difference of the event maximum and pre¬ 
event background spectral indices from -3 to 3, both fast 
and slow injection models, and moderate and weak scat¬ 
tering conditions. The intensity increase timescale was 
obtained from the simulated events by using intensities 
at two consecutive times, with 


Tr = 


In [levih]/Iev{ti)) 


to — t-i 


(9) 


The scaling in Eq. ([^ is obtained by multiplying the 
intensity increase timescale by td/tdo- As the radial 
mean free path A,, is constant in our model, we use the 
Xr and r=l AU in the diffusion timescale. The reference 
value, Tjjo is calculated using Ar = 1 AU and v = c. 

We found that comparing tj with the onset time did 
not result in the expected power law behaviour given by 
Eq. (|^, due to the inability of the diffusion description 
to describe the early time evolution of SEP propagation. 
However, as shown in Pig. the power law behaviour 
is retained when plotting titu/t^q as a function of de¬ 
lay time td, for several orders of magnitude, for different 
particle energies, scattering conditions and injection pro¬ 


files. The delay time is a reminescent of the concept of 
signal speed in modellin g the particle propag ation with 
the Telegrap h equation (Fisk fc Axford ll 19691) . However, 
as shown by Effenberger & Litvinenko (2U14I, the Tele¬ 
graph equation models the initial phase of an SEP event 
poorly, and is not hence addressed further in this stu(W. 

As seen in Fig. the power law suggested by Eq. Q 
is not exact. Thus, to estimate the delw time for using 
it as a correction for VDA, as in Eq. ([^, we fitted the 
data in Fig. to 


td = a 


tiTd 


( 10 ) 


which gave a = 0.61, with standard deviation CTq = 0.15. 
This estimate can be used further in correcting for the 
effects of interplanetary transport on VDA. At larger val¬ 
ues of T/, the diffusive profile is approached and values 
of a up to unity could be used. 

The data can also be fitted as a power law, resulting 
in power law index of 0.56, instead of the form sug gest ed 
by Eq. (§. However, we consider the fit to Eq. ( |l0[ ) a 
better choice, as it is based on the physics of particle 
propagation at the time-asymptotic limit, and as such is 
more likely applicable to other energy ranges and particle 
species. 

We have applied the correction to the observed onset 
times given by Eq. ( 1^ to the six simulated event exam¬ 
ples shown in Fi g. [4 and present the resulting, corre cted 
VDA fits in Fig.|T The delay time td given by Eq. ([lO I 
has been subtracted from the onset times, to, and me 
error limits as defined by Ua have been used in the fit¬ 
ting to obtain error limits for the injection time and the 
path length. As can be seen, the correction improves 
both the injection time and the path length estimates 
considerably. The injection times are correct within the 
error limits for the fast injection profile (left column), 
with the path lengths also better reflecting the Parker 
Spiral length of 1.17 AU. 

When analysing real SEP events, the limited counting 
statistics of the particle detectors cause an uncertainty 
factor for the delay time through the dependence of t/ 
on intensity. We assume that the pre-event background 
intensity can be evaluated from sufficiently long time pe¬ 
riod so that its statistical error is insignificant compared 
to those of intensities at times U and to in Eq. With 
this assumption, the statistical error of t/ can m evalu¬ 
ated as 


Arf 


1 


[ln(/2//i)]^ 


1 

Ad 


1 


where the variables Ni and No represent the number of 
particles detected by the instrument at the two times 
used to calculate t/. Using this approximation, we have 



1 Ni+No 
4iViiV2[ln(/2//i)]' 


( 11 ) 


If the relative error from the background determination 
is significant, it should be included in the error analysis. 

In addition to the Hnite counting statistics, t here is still 
one significant source of uncertainty in Eq. (10). The 


diffusion timescale, td, depends on the scattering mean 
free path of the particles, A. While A depends on the 
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Fig. 6.— Onset times fitted for the same simulated events shown in Fig.j^ with the onset times corrected using Eq. ( |10[ |. The injection 
profiles, parameters, symbols and curves are as in Fig. 
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Fig. 7.— Proton intensities at energies 2.41—90.5 MeV on June 
10, 2000, as observed by the ERNE instrument onboard SOHO 
sp acecraft. The onset times, as det ermined by method described 
in|Huttunen-Heikinmaa et al |||2005|l, are shown with red triangles. 
The intensities in eacP consecutive energy channel are multiplied 
by a factor of 100 to better distinguish onset times. 



Fig. 8.— Event maximum (red curve and squares) and pre-event 
background (blue curve and diamonds) spectra on June 10, 2000, 
with respective power law fits. 


amplitude of interplanetary turbulence, it is difficult to 
estimate, and is typically obtained as a side product of 
SEP transpor t fitting. According to several studies (e.g., 
Palmer|1982 ), the mean free path varies from 0.08 to 0.3 


AU, with some recently analyzed ev ents showing signif - 
icantly longer mean free paths (e.g., Torsti et al. 2004). 
There is no reason to expect the mean tree paths to be 
normally distributed, thus we evaluate its effect on the 
delay time estimation by using extreme mean free path 
values, e.g., 0.1 and 1.0 AU, and obtaining the smallest 
and largest parameter values (taking also the statistical 
errors into consideration). 


3.3. June 10 2000 SEP event 



Fig. 9.— Onset time and velocity dispersion fit correction illus¬ 
tration for June 10th, 2000 event. Observed onset times are repre¬ 
sented by black diamonds, and thin black line is the corresponding 
least-squares linear fit (s = 1.63 di 0.06 AU; to = 16 : 43di4 min). 
Red squares represent the corrected onset times using reference 
mean free path 0.3 AU. See the text for details of the error bars. 
Blue dotted, thick red, and green dash-dotted lines correspond to 
weighted least-squares linear fits for corrected onset times of ref¬ 
erence mean free paths 1.0, 0.3 and 0.1 AU, respectively. These 
three corrected velocity dispersion fits can be incorporated to cor¬ 
rected values of s = 1.03~^S'2ft AU, and to = 16 : 52~’"^^!^. Eight 
- 0.28 ’ ^ -8mm 

minutes need to be added to the times to make them comparable 
to electromagnetic observations. 


In order to study the effects of the delay due to scat¬ 
tering, and the ability of the correction derived above to 
improve the VDA results, we apply the correction to an 
observed SEP event. In Pigs. [TjandJ^ we show the time- 
intensity profiles and spectra m the SEP event of June 10, 
2000, as ob served by ERNE in strument onboard SOHO 
spacecraft (|Torsti et al.||1995|). The event takes place 
during a decay phase of an earlier SEP event. For this 
reason, the pre-event spectrum is very soft compared to 
the event maximum spectral index, as shown in Fig. 
Thus, this event corresponds to the simulated events in 
the bottom row of Fig. 

The SEP event coincides with an M5.2 solar flare 
and a western halo CME. The flare, located at heli¬ 
ographic coordinates N22 W38, started at 16:40 UT 
and reached its maximum at 17:02 UT. The CME was 
launched at 16:45 UT, as given by a linear extrapola¬ 
tion of SOHO/LASCO observations to solar surface in 


the CDAW SOHO/LASCO CME list (Gopalswamy et al., 
2009|. The western location of the buption implies a 


good magnetic connection for the SEPs along the Parker 
spiral to the near-Earth spacecraft, with negligible cross- 
field propagation effects on the first-observed particles. 

We used the VDA method to determine the solar injec¬ 
tion time and path length of the 2.41-90.5 MeV protons 
in this event. Using uncorrected onset times, VDA gave 
for the solar injection 16:43±4 minutes solar time. Tak¬ 
ing into account the 8 minutes light needs to travel from 
the Sun to Earth, the particles are injected 11 minutes 
before the flare maximum time, with the CME estimated 
to be low in the corona. The path length of the parti¬ 
cles obtained from uncorrected onset times is very long, 
s = 1.63±0.06 AU. This is consistent with the simulated 
events (bottom row of Fig. |^: an event with maximum 
spectrum harder than the pre-event spectrum will show 
a very long path length. 

We then applied the correction, as given by Eq. (10), 
to the onset times. We show the effect of the correction 
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in Fig. 1^ where the black diamonds and the black line 
correspond to the uncorrected onset times and the cor¬ 
responding VDA fit. The red squares show the corrected 
onset times, with the associated errorbars, for mean free 
path A = 0.3 AU. We calculated the corrected VDA re¬ 
sults using three different mean free paths, depicted with 
the blue dotted, solid red and green dash-dotted curves 
in Fig. As can be seen in the inset of the Figure, all 
three fits give later time of SEP injection, as compared 
to the VDA fit without time correction. 


Using the errors given by Eq. (11) in the velocity dis¬ 


persion fit, we obtain error limits for each of the fits 
with different A values. Using injection time, path length 
and the fitting error values we find the corrected in¬ 
jection time as -|- 8min, with path length 

s = 1.031 q' 28 Thus, the path length after the correc¬ 
tion is consistent with the expected Parker spiral length. 
The solar injection time, of 16:59 UT corrected for elec¬ 
tromagnetic observations, coincides well with the maxi¬ 
mum time of the X-ray flare. The CME was observed at 
17:08 UT at 2.76 R©. Thus, within the error limits of the 
injection time, both the flare and CME observations are 
consistent with potential energetic particle production. 

4. DISCUSSION 

As this study shows, the results of VDA fitting should 
be used carefully. SEPs scatter in the interplanetary 
medium, and while the first particles related to an event 
may indeed be scatter-free, their intensity is likely too 
low for statistically reliable observation. It should be 
noted that having a very low level of pre-event count¬ 
ing rates, such as in the case of heavy elements, does 
not imply that the first observed SEP event particles are 
scatter-free. The instrument’s detection threshold has a 
similar effect to the observed event onset as the pre-event 
background intensity level. 

It is important to notice that a reasonable path length 
does not imply a good estimate for the solar injection 
time. This can be seen in the fits presented in Fig. 

If the path length is nominal, as in the top panels, ^ 
of the fitted energies will have almost equal delay time 
td, which results in a large error for the injection time. 
If, on the other hand, the path length is long, the de¬ 
lay time is shorter at higher energies (smaller I//3), and 
the VDA-fitted curve converges towards the scatter-free 
VDA pattern (middle panels of Fig. B at the limit of 
1/13 = 0, and the resulting error in tne VDA injection 
time is small. This can take place in particular when the 
spectral indices of the pre-event background and event 
maximum spectra are similar. In this case the ratio of 
maximum and background intensities is independent of 
energy, and the VDA result is simil ar to the metho d 


used by Lintunen & Vainio (2004) and Saiz et al. (2005), 
who determine the onset time as the time when a hxed 
fraction of maximum intensity at the energy channel is 
reached. 

It should be remembered that also factors other than 
the parallel scattering may influence the particle prop¬ 
agation in the interplanetary space. Recent multi¬ 
spacecraft analyses of SEP events suggest that SEP 


(iDresing et al. 20121 iDresing et al. 2014 

Wiedenbeck 

et al. 

2013 Richardson et al. 

sunii^v] 

tide the large 


longitudinal spread of particles may be caused by pro¬ 
cesses low in the corona (see, e.g., Wiedenbeck et aL| 


2013 for discussion), the observed longitudinal depen¬ 
dence of anisotropy suggests that signific ant interplane¬ 
tary cross-field transport is taking place (Dresing et al.j 


2014). While cross-field diffusion has traditionalfy been 


suggested to transport SEPs across the me an field also 


2013 

Marsh et al. 

2013) and propagation along mean- 

dering held lines (Giacalone et al. 

2000 Laitinen et al.| 


Each of the suggested mechanisms would result in dif¬ 
ferent type of interplanetary transport, which may be 
seen in the velocity dispersion pattern. However, with¬ 
out taking the diffusion in the interplanetary medium 
into account in any way, it may be difficult to discern be¬ 
tween the different mechanisms. In this work, we derived 
an estimate for the delay time of SEPs due to the scat¬ 
tering, using a simple model for the onset evolution, and 
SEP transport simulations. We showed that the estimate 
improved the VDA fitting results in the case where par¬ 
ticles propagate along the Parker spiral field. We tested 
the model by analysing a solar event that could be as¬ 
sumed to be magnetically well connected from the Sun 
to Earth and found that the correction brought both the 
injection time and path length to be consistent with the 
expected solar and interplanetary conditions. A compar¬ 
ison of our model results against multi-spacecraft obser¬ 
vations, as well as events modelled with cross-field dif¬ 
fusion, drifts and meandering fieldlines, may bring more 
light into the mechanism behind the efficient spreading 
of SEPs in the inner heliosphere. 


5. CONCLUSIONS 

In this work, we have studied the validity of the veloc¬ 
ity dispersion method in estimating the injection times 
of SEPs by using simulations of energetic protons. We 
find that the typical method of determining the onset 
time as time when the intensity of the SEPs exceeds the 
pre-event background in statistically significant amount 
can lead to significant errors when estimating the solar 
injection time and path length with the VDA method. It 
is important to note that a reasonable path length does 
not imply a good estimate for the injection time. 

We studied the use of the intensity increase timescale, 
Tj to improve the VDA estimate. The gradual increase 
of the SEP intensities in time is caused by the scattering 
of the SEPs in the interplanetary space, and we find a 
relation between r/, the diffusion timescale tjj and the 
delay of the SEPs with respect to the scatter-free propa¬ 
gation time of the particles from Sun to the Earth. Using 
this relation, we showed that the injection time estimate 
given by the VDA can be improved in case of magneti¬ 
cally well-connected SEP events. 

We conclude that the injection times and path lengths 
obtained by using the VDA method should be used with 
care. The interplanetary scattering of particles does de¬ 
lay the arrival of the first-observed SEPs. This delay is 
likely to cause errors in injection analysis based on the 
observed SEP onset times particularly if the pre-event 
particle intensities are high compared to the event max¬ 
imum intensity. While our analysis was performed using 
energetic protons, our results are valid for any particle 
species that are affected by interplanetary scattering, in- 




































































10 


Laitinen et al. 


eluding heavy ions and electrons. 

We acknowledge support from the UK Science 
and Technology Facilities Council (STFC) (grants 
ST/J001341/1 and ST/M00760X/1). The CME catalog 


is generated and maintained at the CDAW Data Center 
by NASA and The Catholic University of America in co¬ 
operation with the Naval Research Laboratory. SOHO is 
a proiect of international cooperation between ESA and 
NASA. 


REFERENCES 


Agueda, N., Lario, D., Vainio, R., Sanahuja, B., Kilpua, E., Sz 
Pohjolainen, S. 2009, A&A, 507, 981 
Aschwanden, M. J. 2012, Space Sci. Rev., 171, 3 
Cane, H. V., Richardson, I. G., & von Rosenvinge, T. T. 2010, 
JGR (Space Physics), 115, A08101 
Dalla, S., Marsh, M. S., Kelly, J., & Laitinen, T. 2013, Journal of 
Geophysical Research (Space Physics), 118, 5979 
Dresing, N., Gomez-Herrero, R., Heber, B., Klassen, A., 

Malandraki, O., Droge, W., Sz Kartavykh, Y. 2014, Astronomy 
Sz Astrophysics, 567, A27 

Dresing, N., Gomez-Herrero, R., Klassen, A., Heber, B., 
Kartavykh, Y., Sz Droge, W. 2012, Sol. Phys., 281, 281 
Droge, W. 2000, ApJ, 537, 1073 

Droge, W. 2003, The Astrophysical Journal, 589, 1027 
Droge, W., Kartavykh, Y. Y., Klecker, B., & Kovaltsov, G. A. 
2010, ApJ, 709, 912 

Effenberger, F., Sz Litvinenko, Y. E. 2014, ApJ, 783, 15 
Fisk, L. A., Sc Axford, W. 1. 1969, Sol. Phys., 7, 486 
Giacalone, J., Jokipii, J. R., Sz Mazur, J. E. 2000, ApJL, 532, L75 
Gopalswamy, N., Xie, H., Yashiro, S., Akiyama, S., Makela, P., 
Usoskin, I. G. 2012, Space Sci. Rev., 171, 23 
Gopalswamy, N., Yashiro, S., Michalek, G., Stenborg, G., 
Vourlidas, A., Freeland, S., Sz Howard, R. 2009, Earth Moon 
and Planets, 104, 295 

He, H.-Q., Qin, G., Sz Zhang, M. 2011, ApJ, 734, 74 
Huttunen-Heikinmaa, K., Valtonen, E., Sz Laitinen, T. 2005, 
ASzA, 442, 673 

Jokipii, J. R. 1966, ApJ, 146, 480 
Kahler, S. 1994, ApJ, 428, 837 

Kallenrode, M., Sz Wibberenz, G. 1990, in International Cosmic 
Ray Conference, Vol. 5, International Cosmic Ray Conference, 

229—h 

Kallenrode, M. B. 1993, Journal of Geophysical Research, 98, 
19037 

Kocharov, L., Vainio, R., Kovaltsov, G. A., Sz Torsti, J. 1998, 

Sol. Phys., 182, 195 


Krucker, S., Sz Lin, R. P. 2000, ApJ, 542, L61 
Laitinen, T., Dalla, S., Sc Marsh, M. S. 2013, ApJ, 773, L29 
Laitinen, T., Huttunen-Heikinmaa, K., Sz Valtonen, E. 2010, 
Twelfth International Solar Wind Conference, 1216, 249 
Laitinen, T., et al. 2000, A&:A, 360, 729 

Lin, R. P., Potter, D. W., Gurnett, D. A., Sz Scarf, F. L. 1981, 
ApJ, 251, 364 

Lintunen, J., & Vainio, R. 2004, AScA, 420, 343 
Marsh, M. S., Dalla, S., Kelly, J., Sz Laitinen, T. 2013, ApJ, 774, 4 
Palmer, 1. D. 1982, Reviews of Geophysics and Space Physics, 20, 
335 

Parker, E. N. 1965, Planet. Space Sci., 13, 9 
Reames, D. V. 1999, Space Sci. Rev., 90, 413 
—. 2009, ApJ, 706, 844 

Reames, D. V., von Rosenvinge, T. T., Sz Lin, R. P. 1985, ApJ, 
292, 716 

Richardson, 1. G., et al. 2014, Sol. Phys., 289, 3059 
Saiz, A., Evenson, P., Ruffolo, D., Sz Bieber, J. W. 2005, ApJ, 
626, 1131 

Torsti, J., Kocharov, L. G., Vainio, R., Anttila, A., Sz Kovaltsov, 
G. A. 1996, Sol. Phys., 166, 135 
Torsti, J., Riihonen, E., Sz Kocharov, L. 2004, ApJ, 600, L83 
Torsti, J., et al. 1995, Sol. Phys., 162, 505 
—. 1998, Geophys. Res. Lett., 25, 2525 

Tylka, A. J., et al. 2003, International Cosmic Ray Conference, 6, 
3305 

Vainio, R., et al. 2013, Journal of Space Weather and Space 
Climate, 3, A12 

Valtonen, E., Kecskemety, K., Kunow, H., Kiraly, P., 

Miiller-Mellin, R., Sz Torsti, J. 2001, J. Geophys. Res., 106, 
10705 

Wiedenbeck, M. E., Mason, G. M., Cohen, C. M. S., Nitta, N. V., 
Gomez-Herrero, R., Sc Haggerty, D. K. 2013, ApJ, 762, 54 
Zhang, M., Qin, G., Sc Rassoul, H. 2009, ApJ, 692, 109 



